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Few-body systems offer a unique challenge to those interested in deep-inelastic scat- 
tering. Using the deuteron as the simplest and most easily solvable nuclear system we 
review the main physics issues in the valence and sea quark regimes. For the former the 
key issue is the understanding of Fermi motion and binding corrections, and our ability to 
extract the ratio F^/Ff as x — 1. The most recent analysis of deuteron data for 
is consistent with perturbative QCD, but inconsistent with the assumptions common to 
all standard parton distribution fits. This should be taken into account in calculating 
event rates at HERA — especially at large x and . The extraction of also presents 
an opportunity to refine the description of nuclear effects in spin-dependent structure 
functions. In the region of small x one must correct for shadowing and meson exchange 
current effects to explore the flavor structure of sea, which is, in turn, deeply related to 
the chiral structure of the nucleon. 



1. INTRODUCTION 

The study of nuclear deep-inelastic scattering, particularly on the deuteron, began in 
the 1970's, primarily as a source of information on the neutron structure function. There 
was an explosion of interest in the subject following the announcement of the nuclear EMC 
effect in 1983 [Q]. For an introduction to the ideas of deep-inelastic scattering and to the 
nuclear EMC effect we refer to recent reviews P,^. While the question of how the structure 
of the nucleon is modified in medium |^ and how that may be investigated using deep- 
inelastic scattering (DIS) [^] is both timely and important, it is not our prime concern here. 
Instead, we concentrate on the deuteron as a source of information about the structure of 
the neutron and as a laboratory for testing our understanding of nuclear corrections such 
as binding, Fermi motion, shadowing and meson exchange currents (MEC). If we cannot 
understand these phenomena in a system as weakly bound as the deuteron, then we shall 
certainly not understand them in any heavier system. 

In section 2 we briefly review the history of binding and Fermi motion corrections, 
leading to a simple one-dimensional convolution of the nucleon momentum distribution 
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with the free nucleon structure function. Recent work on the off-shell corrections to such 
a treatment is then summarized, before we turn to a very interesting physics question, 
namely the validity, or otherwise, of the perturbative QCD (PQCD) predictions for the 
asymptotic d/u ratio in the proton. This issue also seems to have consequences for the 
analysis of the anomaly at large x and recently observed by ZEUS and HI at HERA. 
The extension to the spin dependent structure function of the neutron, g^, is also outlined. 

In section 3 we turn to nuclear corrections at small x, the region dominated by sea 
quarks. The major surprise in this region in recent years was the discovery of a significant 
violation of the Gottfried sum rule 0, which implied a major difference in the number of 
u and d quarks in the proton — a conclusion supported in graphic fashion by recent Drell- 
Yan data from Fermilab experiment E866 [0. The origin of this asymmetry seems to be 
the pion cloud of the nucleon @,|], which is a non-perturbative consequence of dynamical 
chiral symmetry breaking and which cannot be avoided even in the deep-inelastic regime. 



2. THE DEUTERON AT LARGE x 



The early work on DIS from the deuteron was aimed primarily at extracting the neutron 



structure function |]10|-[1^. Amongst the many approaches to this problem we mention the 
light-front treatment |T3| , p^ and the relativistic impulse approximation involving 
the free nucleon structure function at a shifted value of x or ||T7|JT8|1 . A more phe- 
nomeno logical approach, developed by Frankfurt and Strikman |TI5], attempts to derive 
the nuclear correction in the deuteron by extrapolation from higher A as a function of 
the "effective density" of the nucleus. The experimental extraction of is usually made 
using either the phenomenological effective density approach or the older, "pre-EMC" 
theoretical treatments. We shall see that in particular the latter are not satisfactory, and 
that taking into account the main lessons of the EMC effect leads to a surprisingly large 
change in the conclusions concerning the ratio F2 /F2 in the region x — >■ 1. 

The traditional impulse approximation treatment of the deuteron structure function 
assumes a convolution of the free nucleon structure function, ^2^, with the non-relativistic 
momentum distribution /at/d of nucleons in the deuteron, calculated in terms of a non- 
relativistic wave function of the deuteron and its binding energy. Although the binding 
energy is very small, the kinetic energy of the recoiling nucleon plays a significant role 
in forcing the struck nucleon further off-shell than one would usually expect — as first 
pointed out by Dunne and Thomas |T^: 
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where p is the off-shell nucleon momentum, y is the light-cone nucleon momentum fraction, 
and ipD is the deuteron wave function. Note that the "flux-factor", y, was omitted in early 



work, but as noted by many authors |19,20| it is important in practice. 



2.1. Covariant Formulation 

In order to assess the theoretical reliability of the non-relativistic impulse approximation 
one needs to go beyond the usual assumptions made in the convolution approach [^. In 
particular, the ingredients necessary for a covariant, relativistic description are a covariant 
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DNN vertex with one of the nucleons (the spectator to the hard colhsion) on-mass- 
shell, and an off-shell photon-nucleon scattering amplitude ("off-shell nucleon structure 



function") W, the full structure of which was only recently derived in Ref . ||2^ . 

The covariant DNN vertex was written down many years ago by Blankenbecler and 



Cook p3[ and, at least in the rest frame of the deuteron, can be expressed in terms of 
relativistic deuteron wave functions |2^J2^ . The deuteron structure function can then be 
written in terms of the Dirac trace of the product {Aq + 'j'^Ai^)W, where the combination 
[Aq + 'j'^Ai^) includes all of the information on the deuteron vertex and (on-mass-shell) 
spectator nucleon, summed over spin projections and can be calculated in terms of the 
relativistic deuteron wave function. 

The analysis of Ref. showed that the most general form of the operator W (which 
is a 4 X 4 matrix in Dirac space), consistent with the discrete symmetries and gauge 
invariance, which contributes in the Bjorken limit is: 



W = WoI + Wi ^+W2 I 



(2) 



where the Wi are functions of p^, and p ■ q (q is the virtual photon four-momentum). 
Thus, whereas in the free case the nucleon structure function involves the combination: 



Tr[(^ + M)W] ~ MWo + M^Wi + p ■ qW2, 

the deuteron structure function involves: 

Tr[(Ao + l''A,^)W] ~ AoWo + p ■ A,Wi + q ■ A{W2. 

Clearly then, even in the absence of Fermi motion one finds that in general: 
+ F,\ 



(3) 



(4) 



(5) 



unless, as is usually implicitly assumed only one of the Wi is non-zero p2 . 

Having established that, in principle, the structure function of the bound nucleon cannot 
equal the structure function of the free nucleon, the important question is how big the 
difference actually is in practice. To estimate this, one can construct a simple model |22 



of the so-called hand-bag diagram for an off-shell nucleon, in which the A^-quark vertex 
is taken to be either a simple scalar or pseudo-vector, with the parameters adjusted to 
reproduce the free nucleon structure functions It is then possible to manipulate the 
exact deuteron structure function into the form 
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^2 



nD(conv) 



where the convolution component: 



F2''^'-\x,Q') = W[ 

Z Jx 



dy<l)^/niy)Fi'i-,Q' 



(6) 



(7) 



now involves the relativistic momentum distribution of nucleons in the deuteron 
(i.e. (pN/niy) depends on all components of the relativistic deuteron wave function, u 
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|2^). Note that the argument of the free nucleon structure function F^ 
indicates that we need to find a quark with momentum a; in a nucleon with momentum y. 
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Figure 1. F2 / F2 ratio as a function of x for the model of Refs.||2 
accounts for off-shell kinematics, and the on-shell model of Ref.lT 
Ref.p7|. 



131 (solid) which 
(dotted) — from 



Furthermore, the "non-convolution" correction term 5^"^^F.f can be separated into two 
pieces associated with off- mass-shell corrections at the DN and 7*A^ vertices ||27|| . 

The result of the fully off-shell calculation from Ref.[^] is shown in Fig.l (solid curve), 
where the ratio of the total deuteron to nucleon structure functions (F^/F^) is plotted. 
We also show the result of an on-mass-shell calculation from Ref.[|l^] (dotted curve), which 
has been used in many previous analyses of the deuteron data [0,^]- The most striking 
difference between the curves is the fact that the on-shell ratio has a very much smaller 
trough at a; ~ 0.3, and rises faster above unity (at x ~ 0.5) than the off-shell curve, which 
has a deeper trough, at a; ~ 0.6 — 0.7, and rises above unity somewhat later (at x ~ 0.8). 

The behaviour of the full off-shell curve in Fig.l is qualitatively similar to that found 
by Uchiyama and Saito Kaptari and Umnikov and Braun and Tokarev [0, who 
also used off-mass-shell kinematics. However, these authors did not include the (small) 
non-convolution correction term, 5^"^^F^. The on-shell calculation [13], on the other 
hand, was performed in the infinite momentum frame where the nucleons are on-mass- 
shell. The problem with this approach in the past has been the lack of reliable deuteron 
wave functions in the infinite momentum frame, and only recently have the first steps 
been taken in calculating deuteron wave functions on the light-cone ||3^. In practice 
one has usually made use of the non-relativistic, S- and D-state deuteron wave functions 
calculated in the deuteron rest frame, without accounting for the explicit binding effects 
which should show up in the infinite momentum frame in the form of additional Fock 
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components (e.g., A^A^-meson(s) ) in the nuclear wave function. 

Clearly, a smaller D/N ratio at large x, as in Refs.|2^j2^, implies a larger neutron 



structure function in this region. To estimate the size of the effect on the n/p ratio 
requires one to "deconvolute" Eq.(|^) in order to extract F^. To study nuclear effects on 
the neutron structure function arising from different models of the deuteron, one must 
eliminate any effects that may arise from the extraction method itself. Melnitchouk and 
Thomas [pl|] therefore use exactly the same extraction procedure as used in previous 



EMC and SLAC data analyses, namely the smearing (or deconvolution) method 
discussed by Bodek et al. This method involves the direct use of the proton and 

deuteron data, without making any assumption concerning itself. For completeness we 
briefly outline the main ingredients in this method. 

Firstly, one subtracts from the deuteron data, F2 , the additive, off-shell corrections, 
to give the convolution part, p^'^'^°'^''\ Then one smears the proton data, Ff, 
with the nucleon momentum distribution function (pN/oiy) in Eq.(|^ to give = F2 1 ^p- 
The smeared neutron structure function, F2", is then obtained from: 

= Fa""^™""^ - F|. (8) 

Since the smeared neutron structure function is defined as F^ = F^/S'„, we can invert 
this to obtain the structure function of a free neutron: 



5„(f 



D (conv) 
2 



FVSp) . (9) 



In Eq.(|^), the proton smearing factor, Sp, can be computed at each x from the function 
4>N/D{y), and a parameterization of the Ff data ||3^. The neutron structure function may 



then be derived iteratively from Eq.(|^). Taking as a first guess Sn = Sp, the values of F^ 
are smeared by the function (pN/niy), and the results used to obtain a better estimate for 
Sn- The new value for Sn is then used in Eq.(^ to obtain an improved estimate for F2, 
and the procedure repeated until convergence is achieved. 

The results of this procedure are presented in Fig.2, for both the off-shell calculation 
solid) and the on-shell model [r^| (dotted). The increase in the off-shell ratio at 



large x is a direct consequence of the deeper trough in the F^ / F2 ratio in Fig.l. To 
illustrate the role of the non-convolution correction, 5^"^^F2 ■, we have also performed the 
analysis setting this term to zero, and approximating F2 by F2 ^'^°^''\x). The effect 
of this correction (dashed curve in Fig.2) appears minimal. One can therefore attribute 
most of the difference between our analysis and the earlier, "pre-EMC" , on-shell results to 
kinematics — since both calculations involve essentially the same deuteron wave functions. 

2.2. A Test of Perturbative QCD 

The precise mechanism whereby the SU(6), spin-flavor symmetry of the parton dis- 
tributions of the nucleon is broken is a question of fundamental importance in hadronic 
physics. In a world of exact SU(6) symmetry, the wave function of a proton, polarized 
say in the +z direction, has the form 



b T) = 1^ t {ud)s=Q) + \u t {ud)s=i) - ^ 1^ i {ud)s=i) 

{uu)s=i) - ^\di {uu)s=i) , (10) 
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Figure 2. The ratio F2/F2 as a function of x, for the off-shell model (solid), off-shell 
model without the convolution-breaking term (dashed), and the on-shell model (dotted). 
On the right-hand axis we have marked the x ^ 1 limits of the SU(6) symmetric model 
(2/3), and the predictions of the models of Refs. [^, 



;i/4) and |8|,|3l (3/7). 



where the subscript S denotes the total spin of the two-quark component. In this limit, 
the nucleon and A isobar would be degenerate. In deep-inelastic scattering, exact SU(6) 
symmetry would be manifested in equivalent shapes for the valence quark distributions of 
the proton, which would be related simply by uv{x) = 2dv{x) for all x. For the neutron 
to proton structure function ratio this would imply: 



F^ 2 



F 



p 



SU(6) symmetry. (11) 



Of course. Nature does not respect exact SU(6) spin-flavor symmetry. The nucleon 
and A masses are split by some 300 MeV, and empirically the d quark distribution is 
softer than the u. The correlation between the mass splitting in the 56 baryons and the 
large-x behavior of F2/F2 was observed some time ago by Close |^ and Carlitz 



Based on phenomenological ||35| and Regge |36] arguments, the breaking of the symmetry 
in Eq.(^) was argued to arise from a suppression of the "diquark" configurations having 
S = 1 relative to the S = configuration. Such a suppression is, in fact, quite natural 
if one observes that whatever mechanism leads to the observed — A splitting (e.g. 
color-magnetic force, instanton- induced interaction, pion exchange), necessarily acts to 
produce a mass splitting between the possible spin states of the spectator pair, {qq)s, 
with the 5 = 1 state heavier than the S = state by some 200 MeV From Eq. ([T0|) , 
a dominant scalar valence diquark component of the proton suggests that in the x —>■ 1 
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limit F2 is essentially given by a single quark distribution (i.e. the u), in which case: 

— ^ — >0 ; 5 = dominance. (12) 

F2 A u 

This expectation has, in fact, been built into all phenomenological fits to the parton 
distribution data. 

An alternative suggestion, based on perturbative QCD, was originally formulated by 
Farrar and Jackson There it was argued that the exchange of longitudinal gluons, 
which are the only type permitted when the spin projections of the two quarks in {qq)s 
are aligned, would introduce a factor (1 — s)^/^ into the Compton amplitude — in com- 
parison with the exchange of a transverse gluon between quarks with spins anti-aligned. 
In this approach the relevant component of the proton valence wave function at large x 
is that associated with states in which the total "diquark" spin projection^ Sz, is zero. 
Consequently, scattering from a quark polarized in the opposite direction to the proton 
polarization is suppressed by a factor (l — x) relative to the helicity-aligned configuration. 



A similar result is also obtained in the treatment of Brodsky et al. I^Q] (based on 
counting-rules), where the large-x behavior of the parton distribution for a quark polarized 
parallel {AS^ = 1) or antiparallel {AS^ = 0) to the proton helicity is given by: q^^{x) = 
(1 — x)^"~^+^'^^, where n is the minimum number of non-interacting quarks (equal to 
2 for the valence quark distributions). In the x 1 limit these arguments, based on 
PQCD, suggest: 

— ^ ^ ^ ~ '1 Sz = dominance. (13) 

F2 7 u 5 

Note that the d/u ratio does not vanish in this case. Clearly, if one is to understand 
the dynamics of the nucleon's quark distributions at large x, it is imperative that the 
consequences of these models be tested experimentally. 

As explained above, our analysis of the SLAC data points involves no assumption 

whatsoever about F^, the only input is the nucleon momentum density in the deuteron, 
4>N/D{y)- The results of this reanalysis are shown in Fig. 3, at an average value of ~ 12 
GeV^. The data represented by the open circles have been extracted with the on-shell 
deuteron model of Ref.[T^, while the filled circles were obtained using the off-shell model 
of Refs. p2| , p7| . Most importantly, the F2/F2 points obtained with the off-shell method 
appear to approach a value broadly consistent with the Farrar- Jackson prediction of 
3/7, whereas the data previously analyzed in terms of the on-shell formalism produced a 
ratio tending to the lower value of 1/4. 

The d/u ratio, shown in Fig. 4, is obtained by inverting F2/F2 in the valence quark 
dominated region. The points extracted using the off-shell formalism (solid circles) are 
again significantly above those obtained previously with the aid of the on-shell prescrip- 
tion. In particular, they indicate that the d/u ratio may actually approach a finite value 
in the x ^ 1 limit, contrary to the expectation of the model of Refs. [l35| , |36| , in which d/u 
tends to zero. Although it is a priori not clear at which scale these model predictions 
should be valid, for the values of corresponding to the analyzed data the effects of 
evolution are minimal. 

Naturally, it would be preferable to extract F2 at large x without having to deal with 
uncertainties in the nuclear effects. In principle this could be achieved by using neutrino 
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Figure 3. Deconvoluted F^/Ff ratio extracted from the SLAC p and D data p8| , ^ using 
the model of Ref . [P^ , p7[] (sohd circles) and Ref.[|l^ (open circles). 



and antineutrino beams to measure the u and d distributions in the proton separately, 
and reconstructing from these. Unfortunately, as seen in Fig. 4, the neutrino data from 



the CDHS collaboration do not extend out to very large x {x > 0.6), and at present 
cannot discriminate between the different methods of analyzing the electron-deuteron 
data. 

2.3. Relevance of the d/u Ratio to the HERA Anomaly 

The HI and ZEUS experiments at HERA have recently produced a small number 
of events at enormously high which have generated tremendous theoretical interest 
42,43]. For > 10, OOOGeV^ and x > 0.45 the valence parton distributions are cal- 
culated to drop dramatically. The HERA anomaly is essentially the excess of observed 
"neutral current" events (i.e., events of the type e~^p —>■ e~^X) over expectations by roughly 
an order of magnitude. Many exotic explanations of this excess have already been sug- 
gested, indeed, the number of possibilities currently exceeds the number of events. How- 
ever, before the new physics can be worked out one must be sure that the input parton 
distributions used to estimate "background" rates are reliable. 

One glaring problem with the current treatment of the partonic "background" is that all 
of the standard distributions used are constructed to satisfy rf/w— i>Oasx— s>lat low-Q^. 
As we have seen, the recent re- analysis of the deuteron data leads to a d/u ratio which 
appears to be consistent with the prediction of PQCD that d/u— >-l/5a.sx^l. In the 
light of this result it is not only vital to find alternative, more direct measurements oi d/u 
at large x, but those generating standard sets of parton distributions should at the very 
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Figure 4. Extracted d/u ratio (see Fig. 2). Also shown for comparison is the ratio extracted 
from neutrino measurements by the CDHS collaboration 



least present alternative parameter sets consistent with the new analysis of the deuteron 
data. Until parameter sets are constructed which are consistent with d/u ^ l/Sasx— >■ 1, 
at ~ 10 GeV^, one cannot be sure of the reliability of "background" rate estimates 
at the extreme values of and x being probed at HERA. Certainly one cannot hope 
to probe other exciting aspects of hadron structure, such as intrinsic charm and beauty 
T4|-B^, until the d/u issue has been resolved. 



2.4. Spin Dependent Structure Functions 

The spin structure functions of the nucleon, g{^^\ are of tremendous interest at present. 
Experimentally, gi is proportional to the difference of DIS cross sections for ep scattering 
with beam and target helicities aligned and anti-aligned Within the parton model 
it may be written in terms of the parton helicity (loosely "spin") distributions, Ag(x) 



g^], with g^*^^-* the number density of quarks with helicity parallel (anti- 



[q" - + q 
parallel) to the helicity of the target proton: 

1 



(14) 



Intense interest in the spin structure functions began in 1988 when EMC announced a 
large violation of the Ellis- Jaffe sum rule which relates Tp{Q'^) = gi{x,Q'^)dx to 



(3) (S) 

the isovector and octet axial- vector coupling constants, g\ and g)^ . The failure of this 
sum rule, which is not a rigorous consequence of QCD, led to questions about the Bjorken 
sum rule, which relates Tp — r„ to g^^ /Q (modulo QCD radiative corrections [^) and is 
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a strict consequence of QCD. To determine r„ one must measure gi{x), which requires 
a polarized nuclear target such as ^He or D. At present, all neutron data extracted from 
the deuteron are obtained by applying a simple, non-relativistic prescription to correct 
for the D-state component (probability ujd) of the deuteron wave function ||50| : 

g-{x)=(l-^-uJoy\^{x)-g^,{x). (15) 

As we explain below, exactly the same techniques as described in the previous section 
may be used to test the accuracy of Eq. (|T5|) . 

While most interest has been focussed on the issue of sum rules, we stress that the 
shapes of gi{x) and gi{x) contain even more important information. For example, the 
same arguments that led to different conclusions about the behaviour of c?/m as x — >■ 1 also 
give quite different predictions for g^ and as x — 1, namely 1/4 and 3/7, respectively. 
Quite interestingly, while the ratio of the polarized to unpolarized u quark distributions 
is predicted to be the same in the two models: 

Au 

1 ; 5 = or 5*2 = dominance, (16) 

u 

the results for the (i-quark distribution ratio differ even in sign: 

— ; S = dominance, (17) 
a 3 

^1 ; Sz = dominance. (18) 



Using the same techniques described in section |2.1| for the unpolarized case, Mel- 
nitchouk. Filler and Thomas derived the most general, antisymmetric, Dirac tensor 
operator of twist 2 for an off-mass-shell nucleon (see also [0): 

= [p0{fhf,Gp+ 4-i,G,) + 7/375G-,] . (19) 

Once again, one finds that three functions, G,, can be constructed in terms of scalar and 
pseudo-scalar vertices. However, in this case there is a new feature, the function Gp does 
not contribute for a free nucleon, whereas it does contribute in a nucleus. 

The spin-dependent deuteron structure function is given by the trace of G^^'^ with a spin- 



dependent ND amplitude |51], which can be evaluated using the relativistic DNN vertex 
of Ref . p4| . Surprisingly, Fig.5 shows that the ratio of the convolution approximation to 
the fully off-shell calculation (the dashed curve) is even closer to unity in this case than 
in the spin-independent case. 

The comparison between the solid and dotted curves in Fig.5 shows that Eq.(^) is 
reliable at the 2% level for x below 0.7. However, the excellent agreement in this region 
between (|15|) and the full calculation rehes on a knowledge of the deuteron D-state prob- 



ability. As shown in Ref.|^, a change of ood by 2% (e.g. from 4% to 6%) leads to an 
error of order 10% or more in g^. For x > 0.7, on the other hand, the approximation 
( p!5D fails dramatically. This will be extremely important when testing the predictions of 
PQCD for the x ^ 1 behavior of the polarized distributions in Eqs. (|T6|) -(p!^). 

The use of a polarized ^He target reduces the errors arising from the extraction, but 
in that case we have, as yet, no analogue of the relativistic vertex which permitted an 
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Figure 5. Ratio of deuteron and nucleoli structure functions in the full model (solid), 
and with a constant depolarization factor corresponding to ujd = 4.7% (dotted line). 
The dashed curve is the ratio of Qi calculated via convolution to Qi calculated in the 
relativistic model — from Ref. IST . 



estimate of off-shell corrections in the deuteron case. As the binding energy and density of 
the target nucleus increases we expect these corrections to increase as well. In particular, 
while the new structure function Gp contributes only at order {v/c)^ and is negligible for 
the deuteron, one cannot yet discount the possibihty of Gp playing a significant role in 
nuclei heavier than the deuteron. 

3. THE DEUTERON AT SMALL x 

At decreasing values of Bjorken quarks begin to dominate, the impulse approx- 

imation should eventually break down. One expects that coherent, multiple scattering 
effects become important when the characteristic time scale, 1/Mx, of the DIS process 
becomes larger than the typical average distance between nucleon in the nucleus. This 
occurs typically in the region a; < 0.1. Such effects are certainly seen in the low- a; deple- 
tion of the nuclear EMC ratio in heavy nuclei. Any shadowing in the deuteron should 



therefore produce a depletion in the ratio / F2 at small x |53-58|. 

A deviation from unity of F.^ / F.^ at low x is highly relevant for understanding the 
violation of the Gottfried sum rule, which measures the integrated difference between 
the proton and neutron structure functions. The violation, discovered by the New Muon 



Collaboration in 1990 ||59|, was interpreted as a sizable difference between d and u in the 
proton, but relies on an accurate determination of nuclear corrections to F-f. 

An excess of d over u cannot be understood in terms of PQCD — it is a property of 
the "intrinsic sea" of the nucleon. We briefly review the interpretation of this excess, 
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including the widely accepted explanation in terms of the pion cloud of the nucleon, first 
suggested in 1983 P,|60|. As the pion cloud of the nucleon is an unavoidable consequence 
of dynamical symmetry breaking in QCD this experimental discovery offers important 
insight into the realization of non-perturbative QCD in the structure of the nucleon. 

3.1. Nuclear Shadowing 

The rescattering of a virtual photon from both nucleons in the deuteron is usually 
described within the Glauber multiple scattering formalism. Corrections associated with 
a relativistic deuteron wave function amount to a few percent out of a shadowing correction 
to F2 that is only a few percent in total, and hence can be neglected. At small x, nuclear 
binding and Fermi motion corrections can be neglected, and the total deuteron structure 
function written as: 



In modeling the shadowing correction, 5^^^'^'^^ F2 our approach is to take a two-phase 
model, similar to that of Kwiecinski and Badelek |]53|-^. At high virtuality the inter- 
action of the virtual photon with the deuteron is parameterized in terms of diffractive 
scattering through the double and triple Pomeron, as well as scattering from exchanged 
mesons in the deuteron. On the other hand, at low virtuality it is most natural to apply 
a vector meson dominance (VMD) model, in which the virtual photon interacts with the 
nucleons via its hadronic structure, namely the p°, uj and mesons. The latter contribu- 
tion vanishes at sufficiently high Q^, but for < 1 GeV^ it is, in fact, responsible for 
the majority of the dependence on Q^. 

For the diffractive component, Pomeron {IP) exchange between the projectile and tboth 
constituent nucleons models the interaction of partons from different nucleons within the 
deuteron. Assuming factorization of the diffractive cross section, the shadowing correction 
(per nucleon) to the deuteron structure function from iP-exchange is written as a convo- 
lution of the Pomeron structure function, -F/^, with a distribution function //p, describing 
the number density of exchanged Pomerons: 

5(^^)F2^(a;,Q2) = f dy My) Ff{xjp,Q'), (21) 
where, treating the deuteron non-relativistically, |53|-p7|]: 



My) = fd'krSnik'). (22) 

Here, y = x{l + M^/Q"^) is the light-cone momentum fraction carried by the Pomeron 
{Mx is the mass of the diffractive hadronic debris), and xip = x/y is the momentum 
fraction of the Pomeron carried by the struck quark in the Pomeron. The deuteron form 
factor, Spilk."^), is given in terms of the coordinate space wave functions: 

Sn{k') = / drMr)fU\k\r), (23) 



where jo is a spherical Bessel function. Within experimental errors, the factorization 
hypothesis, as well as the y dependence of fip(y) p3H57[], is consistent with the recent 
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HERA data |^ obtained from observations of large rapidity gap events in diffractive 
ep scattering. These data also confirm previous findings that the Pomeron structure 
function contains both a hard and a soft component: Fl^{xip^ Q'^) = F2^^^^'^^\xjp, Q"^) + 
p IP (soft) ^^^^^ Q2y rjj-^g hard component of F^^ is generated from an explicit qq component 
of the Pomeron, and has an xjp dependence given by x/p(l — xjp) in agreement with 



the recent diffractive data The soft part, which is driven at small xjp by the triple- 
Pomeron interaction , has an xip dependence typical of a sea-quark distribution, with 
normalization fixed by the triple-Pomeron coupling constant. 



The dependence of F2 on Q [|^, at large Q\ leads to a weak (logarithmic) Q 
dependence for the shadowing correction 5^^^^F^. The low-Q^ extrapolation of the qq 
component is parameterized by applying a factor Q'^/{Q'^ + Ql), where Ql ~ 0.485 GeV^ 



may be interpreted as the inverse size of partons inside the virtual photon ||63|. For the 
nucleon sea quark densities relevant for j?^^^'^"^^^ -^^e use the recent parametrization of 
Donnachie and Landshoff, which includes a low-Q^ limit consistent with the real photon 
data, in which case the total Pomeron contribution S^^^^F^ — as — > [|63| . 

An adequate description of shadowing at small requires a higher-twist mechanism, 
such as vector meson dominance. VMD is empirically based on the observation that some 
aspects of the interaction of photons with hadronic systems resemble purely hadronic 
interactions. In terms of QCD this is understood in terms of a coupling of the photon 
to a correlated qq pair of low invariant mass, which may be approximated as a virtual 
vector meson. One can then estimate the amount of shadowing in terms of the multiple 
scattering of the vector meson using Glauber theory. The corresponding correction (per 
nucleon) to the nuclear structure function is: 



where 

SCTVD 



2 

' 87r2 



J d'kTSoik') (25) 



is the shadowing correction to the vector meson-nucleus cross section, fy is the photon- 
vector meson coupling strength, and My is the vector meson mass {V = p^,uj,(j)) are 
important at low Q^. The vector meson propagators in Eq.(p3) lead to a strong 



dependence of ^(^^Ff , which peaks at _ 1 GeV^. For ^ and fixed x, 6^^^F^ 
disappears due to the vanishing of the total F^. Furthermore, since this is a higher 
twist effect, shadowing in the VMD model dies off quite rapidly between ~ 1 and 10 
GeV^, so that for > 10 GeV^ it is almost negligible — leaving only the diffractive 
term, S^^^^F^. (Note that at fixed u, for decreasing the ratio F^ / F2 approaches the 
photoproduction limit.) 

3.2. Meson Exchange Currents 

Meson exchange currents (MEG) are familiar in few-body physics. In the context of 
DIS they have been recognized as a potentially important, many-body correction, which 
does indeed scale, since the possible enhancement of the nuclear pion field was proposed 



54i|65| as an explanation of the nuclear EMG effect. In the deuteron case, Kaptari et 
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al. have suggested that MEC may lead to some antishadowing corrections to F,p{x 



The total contribution to the deuteron structure function from meson exchange is: 
5WF2^(x,g2) = J2[dyfMiy)Fi'ix/y,Q'), (26) 



where M = 7T,p,uj,a. The virtual meson structure function, , one approximates by 
the (real) pion structure function, for which data has been obtained through the Drell-Yan 
process. The exchange-meson distribution functions, fhiiy)-, are obtained from the non- 
relativistic reduction of the nucleon-meson interaction given in Refs. [^|66| . In practice 



pion exchange is the dominant process, and gives a positive contribution to 6^^^^F^. The 
exchange of the fictitious a meson (which represents correlated 27i exchange) also gives rise 
to antishadowing for small x. Vector meson exchange cancels some of this antishadowing, 
however the magnitude of this is smaller. In fact, for not too hard meson-nucleon vertices 
all contributions other than that of the pion are essentially negligible. 

3.3. Neutron Structure Function at Small x 

In the region x < 0.1 the magnitude of the (negative) Pomeron/VMD shadowing is 
larger than the (positive) meson-exchange contribution, so that the total S^^^'^'^^F.p = 
^mpD §{v)pD §{M)pD -g negative. On the other hand, at slightly larger x (~ 0.1- 
0.2), there is a very small amount of antishadowing, which is due mainly to the VMD 
and pion-exchange contributions. In the kinematic region covered by NMC, x > 0.004, 
= 4 GeV^ p^, the overall effect on the shape of the neutron structure function is a 
1-2% tncrease in ^27^2"^^°""^^ for x < 0.01, where ^2"^^°""'^) = F^ - Fl 

The presence of shadowing in the deuteron would be confirmed through observation of a 
deviation from unity in the F2 / F2 structure function ratio in the kinematic region where 
Regge theory is expected to be valid. Although the exact value of x below which the proton 
and (free) neutron structure functions become equivalent is not known, it is expected that 
at low enough x, F^ F^, in which case F^^/F^ ^ 1 + S'-''^'"^^F^ /2F^. Subtracting the 
shadowing correction from the deuteron data should then yield — F2 — > as x — > 0. 
While for the lowest NMC data point it may be debatable whether the Regge region is 
reached, the E665 Collaboration |68] has taken data to very low x, x ^ 10^^, which should 



be much nearer the onset of Regge behavior. 

In Fig. 6 we show the ratio F^/Ff extracted from the low x E665 Collaboration data 
as well as the earlier NMC data at larger x which have been corrected for the 



effects of shadowing and MEC |69| (note that the E665 data are not taken at fixed Q^). 
With these corrections, the data are in fact consistent with the hypothesis that F2 = Ff 
in the small x region (say x < 10~^). 

3.4. The Gottfried Sum Rule 

The main interest in the NMC measurement of the neutron structure function at low 
X was an accurate test of the Gottfried sum rule: 

= /',,£M^^ (27) 

JO X 

1 

dx {d{x) ~u{x)). (28) 
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Figure 6. x dependence of the F^f/Ff structure function ratio, compared with the values 

and NMC data 



extracted from the E665 38 



67], and corrected for nuclear shadowing. 



In the naive quark model, where d = u, this gives Sg = 1/3. Ignoring nuclear effects, 
the most recent experimental value obtained by the NMC (over the range 0.004 to 0.8) is 
= 0.2281 ± 0.0065 |T^, indicating a clear violation of SU(2) flavor symmetry in the 
proton's sea. However, because the structure function difference in Sg is weighted by a 
factor 1/x, any small differences between F2 and ^^"■(^°"'^'^) g^^e amplified for x ^ 0. 
Including the nuclear corrections, the overall effect on the experimental value for Sq'- 



Sr 



^exp 
G 



1 ^(shad)^D/ 

dx — 

X 



(29) 



is a reduction of between -0.010 and -0.026, or about 4 and 10% of the measured value 
without the shadowing correction . Therefore a value that reflects the "true" Gottfried 
sum should be around Sg ~ 0.2, which represents some 30% reduction from the naive 
quark model prediction, 5*^ = 1/3. 

3.5. SU(2)j7 Violation and the Chiral Structure of the Nucleon 

The fundamental degrees of freedom in QCD are quarks and gluons, and for a consider- 
able time there was great reluctance to include any other degrees of freedom in modeling 
hadron structure. On the other hand, extensive studies of non-perturbative QCD have 
shown that the chiral symmetry of the QCD Lagrangian is dynamically broken and that 
the resulting, massive constituent quarks must be coupled to pions. As pseudo-Goldstone 
bosons, the latter would be massless in the chiral limit. Most importantly, as emphasised 
by extensive work on chiral perturbation theory, no perturbative treatment of qq creation 
and annihilation can ever generate the non-analytic behaviour in the light quark mass 
for physical quantities (such as M, < >,o't,n) which are generated by these Goldstone 
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bosons. As a consequence, it is now difficult to imagine a realistic quark model which 
does not incorporate at least the pion cloud of the nucleon. 



1^ 




0.3 0.4 



Figure 7. Comparison of the value oi d/u extracted from the recent E866 data [0] with the 
expectations from just the ttN Fock component of the nucleon wave function, for three 
different form factors at the t^N vertex. 



The early development of chiral quark models began in the late 1970's ||60|j7lH 73| and 
it became clear that the pion field had very important practical consequences for the 
low energy properties of hadrons. For example, the coupling to a "bare" nucleon lowers 
its mass by several hundred MeV, and can contribute a large part of the N — A mass 
difference, while the charge form factor of the neutron is a first order effect of its pion 
cloud p]|. 

However, the relevance of the pion cloud for deep inelastic scattering, which had been 
first realized by Feynman and Sullivan [Q, was not explored in depth until the discovery 
of the nuclear EMC effect, when it was suggested that the effect might be caused by 
a nuclear enhancement of the pion field of the bound nucleons [34,65|. At this time it 
was realized that the light mass of the pion would lead to an enhancement of the non- 
strange over the strange sea of the nucleon f^, which was used to put a constraint on the 
pion-nucleon form factor — a limit that has since been explored in detail [[7^. A more 
important consequence of the nucleon's pion cloud, which was also pointed out in Ref. [§], 
was an excess of d over u quarks. In particular, simple Clebsch-Gordan coefficients for 
isospin show that the pion cloud of the proton is in the ratio 2:1 for 7t~^ : 7r°. Since the 
TT"*" contains only a valence d and the 7r° equal amounts of d and u, this component of the 
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pion cloud of the nucleon yields a ratio for d : u of 5:1. 

On the other hand, perturbative QCD inevitably leads to the conclusion that u = d, a. 
result known as SU(2)-flavor, 811(2)^7, symmetry. The latter is a very misleading name 
as there is no rigorous symmetry involved. Indeed, the prediction that the component of 
the sea arising from the long-range piece of the pion cloud of the nucleon satisfies d : u 
= 5:1 is totally consistent with charge independence. Thus the violation of the Gottfried 
sum rule and the subsequent measurement of d — u gives us direct evidence that there 
is a sizable non-perturbative component of the nucleon sea. Various studies of the pion 
cloud of the nucleon since the original NMC measurement have concluded that this is 
indeed the most likely explanation of the observed violation of the sum rule P,|76|- f7^] . 

In Fig. 7 we show the calculated ratio of d/u from the vrA^ Fock component of the 
nucleon wave function as a source of asymmetry, in comparison with the data points 
extracted from the preliminary results [0 from the E866 collaboration on the ratio of 
pD and pp Drell-Yan cross sections. The calculation was performed in the light-cone 
formalism [^,0] with a monopole irNN form factor mass parameter A = 0.7, 1.0 and 1.3 
GeV — from smallest to largest. Clearly the agreement is qualitatively excellent, but a 
more detailed analysis needs to be carried out once all the data have been analysed. Let 
us emphasise once more the importance of these data, which are giving us direct insight 
into the way dynamical chiral symmetry breaking is realized in the nucleon. 

3.6. Other Non-perturbative Components of the Sea 

The experience with the breaking of SU(2)ir symmetry which we have just described 
leads us to take much more seriously the possibility of an intrinsic component of the 



strange and charm quark sea. The former was first discussed by Signal and Thomas |]82 



and has recently been investigated by a number of authors P4,p3[ in the light of new 



neutrino data from CCFR This will clearly be the subject of much more detailed 



investigation in future. With regard to the question of intrinsic charm |4^J46|] , we refer to 
our earlier discussion of the HERA anomaly for some sense of the issues involved. This 
too will clearly be a subject of intense investigation over the next few years. 



4. CONCLUSION 

At this stage it should be clear that the study of deep-inelastic scattering from few- 
body systems, notably the deuteron, provides a rich source of new information on hadron 
structure. It is a fascinating source of challenges for few-body physicists, including such 
issues as relativity, off-shell effects, meson exchange currents and vector meson dominance. 

As the result of a careful analysis of binding effects in deep-inelastic scattering from the 
deuteron, there has been a dramatic change in the experimental value of the d/u ratio in 
the valence region. It seems to be consistent with the perturbative QCD prediction that 
d/u — >■ 1/5 as a; 1, rather than — as assumed in all standard parton distributions 
(and text books). This analysis strongly indicates a need for alternative determinations of 
d/u — e.g. through semi-inclusive measurements — just to be sure. Meanwhile, there is 
an urgent need to work out the quantitative consequences for this finding in the kinematic 
region relevant to the HERA anomaly. 

It was reassuring that the model dependent analysis of the "non-convolution" correc- 
tions in the deuteron gave values of order 1-2%. However, it remains an open question 
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how big they could be in heavier nuclei. In particular, it would be of great interest to 
have a relativistic '^He-N vertex of the Gross type. This is even more important because 
^He is usually used in the determination of g^. 

In the small-x region a careful analysis of recent CERN and Fermilab data in terms 
of VMD, Pomeron and meson exchange corrections suggests that and Ff are equal 
below X ~ 10^'^. A very important discovery in this region has been the violation of 
the Gottfried sum rule. The physics involved in this violation has been made explicit by 
recent data from the E866 Collaboration at Fermilab, which shows the x-dependence of 
the difference between d and u. 

The theoretical origin of the d-u difference seems to be the pion cloud of the nuclcon, 
which in turn is a direct consequence of dynamical chiral symmetry breaking. Further 
study of the shape of this component of the intrinsic sea of the nucleon should throw 
new light on the way chiral symmetry is realized in non-perturbative QCD. From another 
point of view, the verification of this prediction for the light quark sea leads us to be 
bolder in looking for differences between s and s, and c and c. 

To conclude, we have seen that deep inelastic scattering is potentially a gold mine of 
information concerning the structure of hadronic matter. Thus far we have really only 
begun to scratch the surface. 
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